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Foreword 


Tliis  repoil  was  prepared  as  part  of  the  Markov  Models  of  Random  Urinalysis  Sampling 
Procedures  project  (work  unit  92PODD91 1)  and  the  Statistical  Methods  for  Drug  Testing  project 
(program  element  0305889N),  under  the  sponsorship  of  the  Chief  of  Naval  Personnel  (PrRS-63). 
Tire  objectives  of  these  projects  include  effectiveness  of  using  urinalysis  strategies  based  on  time 
since  last  test  to  improve  the  Navy's  drug  deterrence  program  and  development  of  a  unified  set  of 
statistical  methodologies  for  the  analysis  of  drug  testing  programs  and  data.  The  work  described 
hei  e  was  performed  during  FY92  and  FY93.  Related  work  includes  Probability  of  Detection  of 
Drug  Users  by  Random  Urinalysis  in  the  U.S.  A^uv_y  (NPRDC-TN-93-2). 

The  opinions  expressed  in  this  paper  are  those  of  the  authors,  arc  not  official,  and  do  not 
necessarily  reflect  the  views  of  the  Navy  Department. 

Our  thanks  to  Dr.  Jules  Borack  for  improvements  to  this  report. 


MURRAY  W.  ROWE 
Director,  Manpower  Systems  Department 
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Summary 


Background 

Since  aiiy  drug  abuse  impacts  readiness,  he<ilth,  and  safety,  continuing  evaluation  and 
improvement  of  the  Navy's  piogram  is  required.  One  method  of  improving  the  Navy's  program  is  to 
develop  and  analyze  alternative  testing  strategies. 

The  Nuclear  Regulatory  Commission  (NTiC)  proposed  a  urinalysis  testing  strategy  based  on  time 
since  last  test.  That  is,  the  probability  of  a  person  being  tested  depends  on  the  amount  of  time  since 
the  person  was  last  tested.  Southern  California  Edison  (SCE)  has  implemented  a  variation  of  the  NRC 
proposal.  Urinalysis  testing  strategies  based  on  time  since  last  test  are  defined  by  a  higfi  testing  rate 
for  personnel  not  yet  tested  in  a  given  time  period  and  a  low  testing  rate  for  previously  tested 
personnel  with  negative  results  in  a  given  time  period. 

Objective 

The  objective  of  this  woik  is  to  determine  if  urinalysis  strategics  proposed  by  NRC  and 
implemented  by  SCE  could  be  used  to  improve  the  Navy's  drug  screening  program. 


Approach 

The  probability  distribution  of  random  urinalysis  tests  is  modeled  under  a  general  class  of  age- 
test  urinalysis  strategies.  Age-test  is  a  particular  Markov  chain  with  the  probability  of  being  tested 
defined  as  a  function  of  time  since  last  test.  The  NRC  proposal,  tJte  SCE  program,  and  current  Navy 
practice  can  all  be  modeled  as  age-test  Markov  chains.  Vaiious  age- test  .sttategies  are  anal>'zed. 


Results 


'Fhe  age-test  Markov  chain  was  u.sed  to  analyze  five  different  urinalysis  testing  strategies:  two 
NRC  proposed  alternatives,  the  SCE  process,  and  two  models  of  the  Navy's  program  assuming  a  15% 
monthly  testing  rate. 

This  age-test  Markov  chain  has  13  states.  The  process  is  observed  monthly  and  states  1  through 
12  are  defined  by  the  number  of  months  since  an  individual  was  last  tested.  State  13  includes 
individuals  who  were  tested  13  or  more  months  ago. 

Age-test  urinalysis  strategies  can  be  used  to  reduce  the  variance  of  the  number  of  tests  per  person 
per  year,  while  keeping  the  average  number  of  tests  unchanged.  ITiis  allows  a  reduction  in  both  the 
number  of  people  who  are  not  tested  and  the  number  of  people  tested  multiple  times.  The  system  also 
becomes  more  predictable  and,  hence,  potentially  subject  to  gaming  by  dnig  users. 
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The  five  urinalysis  strategics  arc  suniniarizcsl  below. 


Strategy 

Annual 

Testing 

Rate 

(%) 

Constant 

Monthly 

Rate 

Probability  Not  Tested 

Witliin  Witliin  1  Year 

1  Year  Given  Just  Tested 

NRC— A 

103 

No 

.100 

.738 

NRC— B 

300 

Yes 

.032 

.032 

SCE 

130 

No 

.052 

.478 

Navy — A 

180 

Yes 

.142 

.142 

Navy — B 

180 

No 

.081 

.108 

Note.  NRC  -  Nuclear  Regulatory  Commission,  SCE  »  Southern  California  Edison. 


These  strategies  have  widely  varying  annual  testing  rates  (10.^-300%),  and  either  constant  or 
varying  monthly  rates.  The  varying  annual  rates  imply  varying  costs  of  the  programs.  A  preferred 
strategy  should  have  low  annual  testing  rates  to  lower  costs,  a  low  probability  of  not  being  tested 
within  1  year,  and  a  low  probability  of  not  being  tested  within  1  year  given  just  tested.  When  this  last 
probability  is  high,  the  system  is  subject  to  gaming  by  drug  users. 

Conclusions 

Markov  chains  provide  a  framework  for  the  systematic  analysis  of  drug  tesiing  strategies  base»l 
on  time  since  last  test.  The  steady  state  distribution  provides  estimates  of  the  number  of  tests  per 
month  and  the  number  of  f>eople  who  have  not  been  tested  in  the  past  year.  The  distribution  of  the 
number  of  tests  in  a  fixed  time  period  (c.g.,  year),  given  any  initial  state,  can  be  calculated. 
Furthermore,  given  test  or  cost  estimates,  the  relative  merits  of  different  testing  strategies  can  be 
easily  calculated,  hi  general,  age-test  urinalysis  strategies  trade  off  predictability  for  reduced  tail 
area.  Here  we  mean  the  tail  area  of  the  distribution  of  the  number  of  tests  in  a  fixed  time  period. 
Age-test  strategies  provide  fewer  people  not  tested  within  1  year  and  fewer  people  tested 
excessively  during  1  year.  Age-test  strategies  are  also  more  predictable,  have  lower  variance  in  the 
number  of  tests,  and  as  a  result  are  subject  to  t,aining  by  drug  users.  Future  work  will  quantify  these 
trade-offs. 

Both  the  NRC  alternative  A  and  the  SCE  process  have  some  undesirable  properties.These 
strategics  involve  large  differences  in  the  testing  rates  between  people  tested  within  the  past  year 
and  those  who  were  not  tested.  This  implies  that  once  tested  there  is  a  high  probability  of  not  being 
tested  again  within  1  year.  These  probabilities  are  0.74  for  NRC  alternative  A  and  0.48  for  tl'e  SCE 
process.  The  SCE  process  is  such  that  almost  half  the  tests  every  month  are  given  to  people  who 
know  they  will  be  tested.  For  these  reasons  we  do  not  recommend  either  NRC  alternative  A  or  the 
SCE  process. 
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1.0  Introduction 


1.1  Background 


The  U.S.  Navy's  zero  tolerance  drug  policy  has  been  in  effect  since  1981.  Since  then  the  Navy 
has  pursued  an  aggressive  urinalysis  testing  program.  The  objectives  of  this  testing  program  have 
been  to  deter  and  detect  drug  abuse,  as  well  as  provide  data  on  the  prcvaicnce  of  drug  abuse.  All 
officer  and  enlisted  pc-.-sonnel  arc  subject  to  random  urinalysis  testing  on  a  continuing  basis. 
Current  policy  (Chief  of  Naval  Operations  1990),  directs  Navy  commands  to  test  10  to  20%  of  their 
members  each  montli.  The  Navy's  random  urinalysis  patgram  has  been  considered  successful.  The 
proportion  of  service  members  sampled  testing  positive  for  drugs  fell  from  7%  to  1  %  between  1 983 
and  1991.  However,  since  any  drug  abuse  impacts  readiness,  health,  and  safety,  continuing 
evaluation  and  improvement  of  the  Navy's  program  is  required.  One  method  of  improving  the 
Navy's  program  is  to  develop  and  analyze  alternative  testing  strategies. 


The  Nuclear  Regulatory  Commission  (NRC)  (1988)  proposed  a  urinalysis  testing  strategy 
based  on  time  since  last  test.  That  is,  the  probability  of  a  person  being  tested  depends  on  the  amount 
of  time  since  the  person  was  last  tested.  Southern  California  Edison  (SCE)  has  implemented  a 
variation,  (Murray  &  Talley  1988),  of  the  NRC  proposal  at  the  San  Onofre  Nuclear  Generating 
Station.  Urinalysis  testing  strategies  based  on  time  since  last  test  are  defined  by  a  high  testing  rate 
for  personnel  not  yet  tested  in  a  given  time  period  and  a  low  testing  rate  for  previously  tested 
personnel  with  negative  results  in  a  given  time  period.  This  strategy  may  provide  balance  among 
program  objectives,  includmg  detection,  deterrence,  a  high  probability  of  tesung  some  fixed 
number  of  times,  a  low  probability  of  testing  more  than  some  fixed  number  of  times,  cost 
effectiveness,  ease  of  administration,  and  avoidance  of  discrimination.  Tlic  NEC's  (Nuclear 
Regulation  Commission,  1989)  adopted  rules  and  regulations  for  urinalysis  do  not  require  a  time 
since  last  test  strategy.  However,  SCE  continues  to  use  their  variation  of  this  strategy  with  NRC 
approval. 


The  advantage  of  a  time  since  last  test  strategy  is  shown  in  the  following  example.  Assume  die 
Navy's  objectives  at  each  command  are:  90%  of  personnel  should  be  tested  at  least  once  each  year 
and  personnel  already  tested  within  1  ycai'  should  be  tested  at  a  rate  of  2.5%  per  month.  A  simple 
random  sampling  (with  replacement)  strategy  that  meets  these  objectives  wou.ld  require  an  average 
of  2.1  tests'  per  person  per  year.  However,  a  strategy  based  on  time  since  last  test  could  meet  these 
objectives  with  an  average  of  1.0  tests  per  person  per  year.  An  annual  savings  of  5.94  million 
(based  on  a  decrease  of  1 .1  tests  per  person  for  600,()00  fieople  at  $9  pier  test)  could  be  realized  by 
using  an  age-test  strategy  in  this  example. 


Markov  chains  provide  a  framework  for  the  systematic  analysis  of  drug  testing  strategies  based 
on  time  since  last  test.  The  current  Navy  simple  random  sampling  strategy  is  included  within  this 
framework.  A  related  use  of  Markov  chains,  modeling  classes  of  drug  users,  is  given  in  Evanovich 
(1985).  Previous  work  (Thompson  &  Boyle,  1992),  includes  models  of  detection  of  drug  users. 


’Based  on  a  monthly  rale  t,  die  probability  of  being  tested  at  least  once  within  12  months  is  1  -  (1  -  i)’^  ■  0.9. 
Therefore  t  -  0. 1 75  and  the  average  amiual  number  of  tests  is  1  2t  -  2. 1 . 
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1.2  Objective 

The  objective  of  this  work  is  to  determine  if  urinulysis  stratejries  proposed  by  NRC  and 
impieniciitcd  by  SCE  could  be  used  to  improve  the  Navy's  drug  scix'cning  program. 

1.3  Approach 

ITic  probability  distribution  of  random  urinalysis  tests  is  modeled  under  u  general  class  of  aj^e- 
urinalysis  strategics.  Age-test  is  a  particular  Markov  chain  with  the  proba'.jility  of  being  tested 
defined  as  a  function  of  time  since  last  test.  The  NRC  proposal,  the  SCE  program,  and  current  Navy 
practice  can  all  be  nuxleled  as  age-test  Markov  chains  witfi  13  states,  -n  this  case,  states  1  through 
1 2  coirespond  to  1  through  1 2  months  since  last  test  and  state  1 3  is  defined  as  over  1 2  months  since 
last  test.  The  probability  of  being  tested,  given  the  current  siate,  defines  an  agr-icst  urinalysh 
strategy.  Various  age-test  strategics  arc  analyzed. 


2.0  Markov  Chains 


This  section  briefly  develops  the  theory  and  notation  that  will  be  used  in  the  remainder  of  the 
report.  Chapters  1  and  2  of  Hcxrl,  Port,  &  Stone  (1972)  arc  the  primary  source  for  this  materiai. 


Consider  a  system  that  can  be  in  any  one  of  a  finite  number  of  states.  This  set  of  states  is 
denoted  Land  is  called  the  state  space  of  the  system.  The  system  is  observed  at  discrete  points  in 
time  .n  =  0,  1,2, 

random  variable;  that  is,  the  system  is  not  deterministic 


and  X  denotes  the  ®tate  of  the  system  at  time  ft. 


mir  nnmfWf»c  V  it  ti 
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Systems  with  the  property  that  only  the  present  state  influences  the  future  are  called  Markov 
chains.  In  such  systems,  knowledge  of  the  path  taken  to  reach  the  present  state  cannot  help  predict 
the  future.The  Markov  chain  is  a  simple  generalization  of  systems  of  independent  random 
variables. 


The  Markov  property  is  defined  by 

4  1  ~  +  1  I  ^0  ~  "'O’  •  •  •»  4  1  ~  -^n  4i  I  ~ 


The  conditional  probabilities  P(X„  |  =  x„)  are  called  transition  probabilities.  When 

the  transition  probabilities  arc  independent  of  time  n,  they  arc  called  stationary  transition 
probabilities.  The  model  developed  in  this  paper  has  random  variables  that  satisfy  the  Markov 
property  and  have  stationary  transition  probabilities.  We  represent  the  transition  probabilities  by 
the  matrix  P,  where 

Pf.x,  y)  -  -  y  I  Xi  =  x). 


The  product  of  P  with  itself  n  times  yields  the  /i-stcp  transition  matrix  w'hcrc 

l^(x.  y)  =  PfX*,  „  =  >  I X*  =  jfk 
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Given  a  stationaiy  transition  niati  ix  P  and  an  initial  distribution  Fly  all  probabilities  assixiiatcd 
with  thr  chain  arc  uniquely  determined.  In  piu'ticular.  the  distribution  ol  is 

n„-noi»". 

Hitting  times  arc  random  variables  which  {.day  an  important  role  in  the  theory  of  Maikov 
chains  For  y  eL,  the  hitting  time  iy  is  defined  as  the  first  positive  time  the  ehain  resides  at  state 
y,  i.c.. 


Ty  min  (n  >  1  :  A'„  «  y). 


Define 

p„  -  <  oo)  -  PiTy<  oo  u„  -  A) 

as  the  probability  that  a  Markov  chain  staning  at  .v  will  hit  >  in  finite  time.  A  state  y  is  transient  if 
py.y  <  1,  or  if  starting  at  y  there  is  some  positive  probability  that  tire  chain  w'ill  never  return  to  y.  A 
state  is  recurrent  if  Pyy-  1 .  Fuiehcrmotc,  ail  rccuneni  states  in  a  finite  state  Markov  chain  have  the 
property  that  the  mean  return  time  is  positive  and  finite,  where 

»iy-£y(7;  j-£(7^  I  Ao-y). 

Such  states  are  called  positive  rccuneni.  A  fundamental  theorem  of  finite  Maikov  chains  states 
that  the  class  of  positive  recurrent  states,  L^,  is  nonempty  and  is  partitioned  into  closed  irreducible 
subclasses.  A  class  of  slates  C  is  closed  irreducible  if,  once  in  C  a  chaiit  cannot  leave  C  and  all  states 
in  C  lead  to  each  other  w  ith  positive  probability. 

Associated  with  certain  Markov  chains  arc  special  distributions  w'hich  satisfy 

^  n  (a  J  P(.v,  y)  -  n(y),  ye  L . 

A' 


Such  a  distribution  H  is  called  stationary  and  a  chain  with  a  stationary  initial  cxistribulion  will 
have  tlie  property  that 


£(X„-y)-n(y),y  e  L 


for  all  /I  >  0.  When  a  stationary  distribution  H  satisfies 

lim  P"(A.y)  -  n(y).  A,  y  e  L 
«— > » 


then  n  is  called  a  steady  state  distribution.  This  is  a  strong  condition  and  implies  that  FI  is  the 
unique  stationary  distribution  and 

lim  £(A'„  =  y)  “  lim  Hq  W  P”fA.  y)  H  (y),  a,  y  €  L 

fl  — ♦  -30  )  CO 
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lor  an  arbitiaiy  initial  distribution  0,,.  liquation  1  means  that,  rcgaidlcss  ot  the  chain’s  initial 
distribution,  lor  lai;gL  n  the  distribution  ot  A'„  approximates  the  steady  state  distribution.  It  is 
impoilam  to  knov;  when  stationary  distributions  exist  and  when  a  stationaiy  distribution  is  steady 
state.  Theoiein  5  (lUvl  ei  al.,  1972,  p  64)  slates  tliat  an  inedueible  positive  recuirent  Maikov  chain 
has  a  unique  staliop;uy  distribution  given  by 


n(yl  - 


y  GL 


Thcoi'cm  7  (Hoel  ct  al.,  1972,  p.  73)  states  that  when  such  a  chain  is  aperiodic  the  stationary 
distribution  dclincd  by  eqiiatitMi  2  is  steady  state.  A  sufficient  condition  for  apcriiKlicity  is  that 
P(.v, .»)  >  0  for  at  least  one  a  G  Li¬ 
lt  is  al.sonccessiuy'  to  detine  the  concept  of  occup  ition  times  and  their  associated  probabilities. 
Define  the  (K’cupation  time 


as  the  number  of  times  Uic  chain  is  in  state  y  in  n  time  pcriixJs,  where  1  ^  is  the  indicator  lunction 


Xoifzfy. 


The  occupation  time  piiibabilitics  are 

Cv)  =  m  I  Xq  -  x)  (4) 

for  m  -  0, 1 . n.  Ihi-s  is  the  probability  of  ni  visits  ii)  state  y  in  n  time  pcritxls  given  that  Lbe  chain 

started  in  state  a. 

3.0  Age-Tt*st  Model 

Using  the  previously  described  notation  and  theory,  we  now  define  a  class  of  age-test  Markov 
chains  where  the  states  arc  defined  by  time  since  la.st  test,  llic  transition  matrix  is 


9,  «  0 
0  0 
0  0 


0  0  0 
0  0  0 


\ 


A  similai'  mixlcl,  called  agc-rcplaccmcnt,  where  *i  -  1  is  presented  in  Taylor  and 
Karlin  (19S4).  Later  we  will  be  primarily  interested  in  the  ease  d  ~  \2  and  testing  is  conducted 
monthly.  In  general,  an  individual  is  in  state  iif  lact  tested  itiinc  periods  ago.  An  individual  resides 
in  state  d  +  1  if  last  tested  d  +  1  or  more  time  f>criixls  ago.  For  each  state  an  individual  is  either  tested 
wi»h  probability/!,,  in  which  ease  the  individual  moves  to  state  1,  or  the  individual  ages  1  time 
periixl  with  probability  -  1  -  /?,  .Hence,  the  chain  is  named  age-test.  For  example,  the  Navy 
program  with  testing  unrelated  to  time  since  last  test  and  a  15%  monthly  testing  rate  would  have 
all  the  /j,’s  equal  to  0. 1 5. 

When  the  restrictions  0  <r  1  and  0  <  pj*  ,  <  1  are  Splaccd  on  the  transition  matrix  in 

equation  5,  all  states  lead  to  al!  other  states  with  positive  prob.ability.  Hence,  the  chain  constitutes 
a  single  irreducible  positive  recurrent  class  and  there  is  a  unique  stationary  distribution.  Let  FI 
-  (T:,,  Ttj, .  itrf  ^i).  Solving  n  -  nP  with  Z  Ttj;.  -  1  will  give  the  stationary  distribution  FL  Writing 
out  the  equations  we  have  the  following 


/J|7t,  PjTTi 

+  ...  + 

Pd^d 

+  Pd*\^d+\ 

- 

1 

• 

^2^2 

♦  1 

(id^d 

^  Qd*  1  1 

1 

- 

7t|  4  ^2 

+  ...  + 

+  ^d*\ 

Solving  in  terms  of  n,  we  obtain 

•i 

712 

- 

9iJ^i 

- 

92^  '  ^2^1  ^1 

% 

- 

- 

'yd -  1  “?d -  \9d-i-  ■ 

^d*  1 

- 

Pd+\)^J~  ■ 

-  •9i  Va.i)^i- 

Since  -  I  then 

7t.  -  1/(1 +«i+^i92  +  9i^2<?3+  -  +(9i92--9a-i)  +  (9i<?2 


Since  p  >  0,  the  chain  is  aperiodic  and  FI  is  steady  stale  with,  according  to  equation  2,  )i,  -  Vm, . 
This  is  the  reciprocal  of  the  mean  return  time  to  state  i . 
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In  Section  4  wc  calculate  the  (x:cupation  time  probabilities  for  state  1  using  the  age-test  nuxlcl. 
The  occupation  time  probabilities  for  1 2  time  periods  provide  the  distiibution  of  the  annual  number 
of  tests.  A  visit  to  state  1  is  equivalent  to  being  tested.  Observe  that  from  any  initial  state  there  are 
2'^  possible  paths  over  12  periods.  "iTie  probability  of  each  path  is  calculated  and  the  results  are 
aggregated  to  get  the  conditional  probabilities  in  equation  4.  The  unconditional  distribution  of  the 
number  of  tests  in  12  time  pericxls,  (/V,2(l)),  with  initial  distribution  equation  6  is  determined  by 

13 

PrjV,2  ( 1 )  -  m)  -53  P,  (A'l  2  ( U  »  m). 

1-  1 

Wc  also  observe  the  average  annual  number  of  tests  is 

f  \ 

£(A^, 2  (!))-£■ 

Vn-1  j 

12 

in  =  1 
=  12  rt, 


since  FI  is  stationary. 

Finally,  note  the  probability  of  not  ing  tested  within  1  year  is 

W2(l)-0)-n,v  (7) 

Specifically,  tlie  event  ,V|2  ( 1 )  -  0  is  equivalent  to  the  chain  being  in  state  1 3  at  time  1 2.  Because 
the  initial  distribution  is  stationary,  all  distributions  are  stationary  and  P  (X,2  -  1 3)  -  See  Sellke 
( 1 984)  for  further  details. 

4.0  Applications 

The  age-test  Markov  chain  was  used  to  analyze  five  different  urinalysis  testing  strategies:  two 
NRC  proposed  alternatives,  the  SCE  process,  aiid  two  models  of  the  Navy's  program  assuming  a 
15%  monthly  testing  rate. 

This  age-test  Markov  chain  has  13  states.  The  process  is  observed  monthly  and  states  1  through 
12  are  defined  by  the  number  of  months  since  an  individual  was  last  tested.  State  13  includes 
individuals  who  were  tested  13  or  more  months  ago.  Age-test  urinalysis  strategies  can  be  used  to 
reduce  the  vaiiance  of  the  number  of  tests  per  person  per  year,  while  keeping  the  average  number 
of  tests  unchanged.  This  allows  a  reduction  in  both  the  number  of  people  who  are  not  tested  and 
the  number  of  people  tested  multiple  times.  The  system  also  becomes  more  predictable  and,  hence, 
potentially  subject  to  ganting  by  drug  users.  Tlie  specific  examples  follow. 
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4.1  Nuclear  Regulatory  Commission  Proposals  (NRC) 

Tvvo  urir'alysis  testing  alternatives  proposed  by  the  NRC  (1988)  are  analyzed  in  this  section. 
Alternative  A  required  that  at  least  90%  of  the  individuals  are  tested  each  year  and  that  testing  rates  for 
individuals  already  tested  with  negative  results  be  at  least  2  per  month.  The  age-test  model  with 

p,  >0.025, /-I,...,  12 

71,3  <0.1  (8) 

satisfies  these  requirements.  Minimizing'^  tr,,  the  average  number  of  tests  p)er  month  per  person, 
subject  to  equation  8  yields  a  solution  with  p,  - P2  -  •  •  •  “Pi2  “  0.025  and p,3  -  0.6338.  See  the 
Appendix  for  a  detailed  analysis.  Results  are  summarized  in  Table  1  and  Figure  1 .  The  last  column 
in  the  table  and  the  figure  were  calculated  using  equation  4  and  enumerating  all  possible  transitions 
over  12  months.  The  advantages  of  this  alternative  include:  (1)  the  average  number  of  tests  per 
person  per  year  is  1.03,  (2)  90%  of  the  people  are  tested  at  least  once  per  year,  and  (3)  12%  of  the 
people  are  tested  more  than  once.  The  major  disadvantage  is  that  once  tested,  people  have  a  74% 
chance  of  not  being  tested  for  the  next  12  months. 

Table  I 


Probabilities  From  Age-Test  Model  of  NRC  Alternative  A 


State 

Prol>aKilitx 

1 

Tested 
This  Month 

Steady 

State 

Not  Tested 

Next  12  Months 

1 

.025 

.086 

.738 

2 

.025 

.084 

111 

3 

.025 

.082 

.104 

4 

.025 

.080 

.039 

5 

.025 

.078 

.015 

6 

.025 

.076 

.006 

7 

.025 

.074 

.002 

8 

.025 

.072 

.001 

9 

.025 

.070 

.000 

10 

.025 

.068 

.000 

11 

.025 

.067 

.000 

12 

.025 

.065 

.000 

13 

.634 

.100 

.000 

Note.  NRC  «  Nuclear  Regulatory  Commission. 


^This  and  .subsequent  optimization  problems  were  solved  using  the  Microsoft  Excel  Solver,  Microsoft 
Corporation  (1991),  on  an  IBM  compatible  personal  computei. 
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lisit  NRC  Nucleai  Regulation  Commission. 

Figure  1.  Steady  state  distribution  of  number  of  tests  within 
12  months  for  age-test  model  of  NRC  alternative  A. 

Alternative  B  required  hat  tests  must  be  administered  throughout  the  year  at  an  annual  rate 
equivalent  to  300%  of  the  population.  The  age-test  model  with 

Pi  >  0,  /  “  1, . .  .,13 

n,  -  300% /12  months -0.25  (9) 

meets  tltese  requirements.  Letting  Pi  -  P2  “  •  •  -  •  Pi3  “  satisfies  equation  9.  Results  are 
summarized  in  Table  2  and  Figure  2.  The  advantages  of  this  alternative  include:  97%  of  the  people 
are  tested  at  least  once  and  the  cunent  state  provides  no  information  on  the  future  of  the  pixxess. 
Any  strategy  with  all  the  p,  equal  has  the  latter  advantage.  When  the  p,’s  are  equal,  the  distribution 
of  the  number  of  tests  in  I  year  is  simply  the  binomial.  A  disadvantage  of  alternative  B  is  the  high 
cost  of  three  tests  per  person  per  year. 
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Probability 


Table  2 


Probabilities  From  Age-Test  Model  of  NRC  Alternative  B 


State 

Tested 

This  Month 

Probability 

Steady 

State 

Not  tested 

Next  12  Months 

1 

.25 

.250 

.032 

2 

.25 

.188 

.032 

3 

.25 

.141 

.032 

4 

.25 

.105 

.032 

5 

.25 

.079 

.032 

6 

.25 

.059 

.032 

7 

.25 

.044 

.032 

8 

.25 

.033 

.032 

9 

.25 

.025 

.032 

10 

.25 

.019 

.032 

11 

.25 

.014 

.032 

12 

.25 

.011 

.032 

13 

.25 

.032 

.032 

Nfltt.  NRC  -  >iuclear  Regulatory  Commission. 


Note.  NRC  -  Nuclear  Regulatory  Commission. 


Figure  2.  Steady  state  distribution  of  number  of  tests  within 
12  months  for  age-test  model  of  NRC  alternative  B. 
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4.2  Southern  California  Edison  (SCE) 


SCE  has  implemented  a  composite  random  sampling  (Murray  &  Talley,  19H8),  approach  to 
urinalysis.  Their  approach  is  based  on  a  sampling  scheme  that  is  part  sampling  with  replacement 
and  part  sampling  without  replacement.  The  entire  population  is  sampled  at  a  specified  rate  with 
replacement.  People  who  have  not  been  sampled  within  the  past  year  are  sampled  at  another 
specified  rate  witliout  replacement,  SCE  states  a  5%  annual  chance  of  not  being  tested  and  a  130% 
average  annual  testing  rate.  The  process  in  use  at  SCE  can  be  modeled  as  an  age-test  Markov  chain. 
An  age-test  model  with 

Tti  <  130%/ 12  months  -  0.1083 
n,3  <  0.05(10) 

Pi  “  P2^  -  ~P\2 


meets  the  specifications  of  the  SCE  sampling  scheme.  There  is  no  feasible  solution  to  equation  10. 
The  following  two  optimizations  were  performed  to  find  an  age-test  strategy  as  close  to  meeting 
equation  10  as  possible.  Recall  that  jt,,  the  steady  state  ptobability  for  state  1,  represents  the 
average  monthly  testing  rate,  and  Jt;3,  the  steady  state  probability  for  state  13,  represents  the 
average  annual  not  tested  rate.  Minimizing  Ttj  subject  to  ^e  restrictions  on  />,  and  n,,  in  equation 
10  yields  “  P2  ”  ■  •  ■  “  Pu  *"  0.0645,  -  1.0  and  nj  -  0.11 13.  Minimizing  it,,  subject  to  the 

re.suictionsonp.andtti  in  equation  lOyieldsp,  -P2” •  •  •  “Pi2“  0.0596.p,-,-  1.0 and rt,3- 0.05 18. 
The  solutions  from  both  of  these  optimizations  are  close  to  tlie  results  given  by  SCE  (Murray  & 


inOON 


IvCSUltS  frOui  the  latter  OptiuiiZatiOn  SUiiuiidTiZCu  iii  Tabic  3  mm  riguic  /\ 
problem  with  this  process  is  that  every  month  almost  one  half  (0.052/0.108  -  0.48)  of  the  tests  are 
given  to  people  who  are  in  state  13  and,  therefore,  know  they  arc  being  tested  that  month. 


Table  3 


Probabilities  From  Age-Test  Model  of  SCE’s  Process 


Probability 

State 

Tested 

This  Month 

Steady 

State 

Not  Tested 

Next  12  Months 

1 

.06 

.108 

.478 

2 

.06 

.102 

0 

3 

O'" 

.096 

0 

4 

.Ot 

.•)90 

0 

5 

.06 

.085 

0 

6 

.06 

.080 

0 

7 

Oo 

.075 

0 

8 

.06 

.070 

0 

9 

.06 

.066 

0 

10 

.06 

.062 

0 

11 

.06 

.059 

0 

12 

.06 

.055 

0 

13 

1 

.052 

0 

Note.  SCE  -  Southetn  Califomu  Ediica. 
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1.0 


Note.  SCE  »  Southern  Calitbmia  Edison. 


FIpiiiv  3,  idst?  distHli'iiticn  of  nuinb^r  of  tosts  wliHln 

12  months  for  age-test  model  of  SCE’s  process. 

4.3  United  States  Navy  (U.S.  Navy) 

U.S.  Navy  policy.  Chief  of  Naval  Operations  (1990),  directs  commands  to  test  10  to  20%  of 
their  personnel  each  month.  Age-test  mt>dcls  with 

0.10  <  p,  <  0.20,  tl,...,  13  (11) 

meet  this  requirement^.  TWo  different  strategies  consistent  with  equation  11  are  presented  here.  Lt 
the  first  strategy  (Navy-A),  let  pj  «...  «  p,3-  0.15.  Choosing  tlie  p,  all  equal  is  consistent  with 
current  Navy  practice.  In  this  model  testing  is  independent  of  time  since  last  test.  The  value  0.15 
was  chosen  because  it  is  the  midpoint  between  0.10  arid  0.20.  Results  are  summarize  i  in  Table  4 
and  Figure  4. 


^Strictly  speaking  0. 1 5  <  « j  <  0.20  is  sufficient  to  meet  this  requirement. 
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Probability 


Table  4 


Probabilities  From  Agw-Test  Model  of  Navy  Program 
At  15%  Monthly  Testing  Rate  (Navy-A) 


State 

Probability 

Tested 
This  Month 

Steady 

State 

Not  Tested 

Next  12  Months 

1 

.15 

.150 

.142 

2 

.15 

.127 

.142 

3 

.15 

.108 

.142 

4 

.15 

.092 

.142 

5 

.15 

.078 

.142 

6 

.15 

.067 

.142 

7 

.15 

.057 

.142 

8 

.15 

.018 

.14? 

9 

.15 

.041 

.142 

10 

.15 

.035 

.142 

11 

.15 

.030 

.142 

12 

.15 

.025 

.142 

13 

.15 

.142 

.142 

Figure  4.  Steady  state  distribution  of  number  of  tests  within  12  months  for  age-test 
modd  of  Navy  program  at  15%  monthly  testing  rate  (Navy-A). 
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Recall  n^,  the  steady  state  pwbability  of  being  in  state  1,  is  the  average  monthly  testing  rate.  A 
second  scategy  (Navy-B)  with  tt.  equal  to  Navy-A  is  presented  for  comparison  purposes,  'fhe 
followiiig  minimization  problem  was  solved;  minimize  Tt,;.  subject  to  Tii-  0.15  an'*  equation  11. 
Note  Tt,,  is  the  probability  of  not  being  tested  within  1  year.  Therefore,  this  strategy  has  the  same 
average  monthly  testing  rate  and  minimizes  the  probability  of  uOt  being  tested  within  1  year.  The 
solution  is  {p,  -pi  "  0.1,p4-0.1 145,p5  «...  -  p,,  -  0.2).  Results  are  summarized  in  Table 
5  and  Figure  5.  In  this  case  jr,,  -  .81  as  compared  to  7t,j  -  .142  for  Navy-A.  Therefoi-e,  with  an 
equal  number  of  montlily  tests  and  keeping  the  monthly  testing  rate  between  10  and  20%,  strategy 
Navy-B  reduces  the  number  of  people  not  tested  in  a  given  year  by  43%  ([.142  -  .081]/  .142  x  1(X)). 


Table  5 

Probabilities  From  Age-Test  Model  of  Navy  Program  at 
180%  Annual  Testing  Rate 


State 

Probability 

Tested 

This  Month 

Steady 

State 

Not  Tested 

Next  12  Months 

1 

.10 

.150 

.108 

1  i\ 

.  lU 

t  ^  iT 

.iDj> 

3 

.10 

.122 

.086 

4 

.11 

.109 

.076 

5 

.20 

.091 

.069 

6 

.20 

.077 

.069 

7 

.20 

.062 

.069 

8 

.20 

.050 

.069 

9 

.20 

.040 

.069 

10 

.20 

.032 

.069 

11 

.20 

.025 

.069 

12 

.20 

.020 

.069 

13 

.20 

.081 

.069 

13 
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Number  of  Tests 

Figure  5.  Steady  state  distribution  of  number  of  tests  within  12  months  for  age-test 
mode!  of  Navy  program  at  anniai!  testing  rate  (Navy-B). 

The  five  urinalysis  strategies  presented  in  this  section  are  summarized  in  Table  6.  These 
strategies  have  videly  varying  annual  testing  rates  (103  to  300%).  This,  of  course,  implies  widely 
varying  costs  of  the  programs.  A  pieferred  strategy  should  have  low  annual  testing  rates  to  lower 
costs,  a  low  probability  of  not  being  tested  within  1  year,  and  a  low  probability  of  not  being  tested 
within  1  year  given  just  tested.  When  this  last  probability  is  high,  the  system  is  subject  to  gaming 
by  drug  users. 

Table  6 

Sununary  of  Age-Test  Markov  Chain  Analysis  of  Five  Urinalysis  Strategies 


Strategy 

Annual 

Probability  Not  Tested 

Testing 
Rate  (%) 

Within 

1  Year 

Within  1  Year 
Given  Just  Tbsted 

NRC-A 

103 

.100 

.738 

NRC-B 

300 

.032 

.032 

SCE 

130 

.052 

.478 

Navy-A 

180 

.142 

.142 

Navy-B 

180 

.081 

.108 
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5.0  Conclusions 


Markov  chains  provide  a  framework  for  tlic  systematic  analysis  of  drug  testing  strategies  based 
on  time  since  last  test.  Under  mild  conditions  on  the  transition  probabilities,  these  Markov  chains 
converge  to  steady  state.  The  steady  state  distribution  provides  estimates  of  the  number  of  tests  per 
month  and  the  number  of  people  who  have  not  been  tc.stcd  in  the  past  ycai-.  Since  this  steady  state 
solution  can  be  expressed  in  closed  form,  optimization  problems  involving  steady  state  can  be 
formulated  and  solvcd.Tlic  distribution  of  the  number  of  tests  in  a  fixed  time  period  (e.g.,  year), 
given  any  initial  state,  can  be  calculatcd.Furthcrmore.  given  test  or  cost  estimates,  the  relative 
merits  of  different  testing  strategies  can  be  easily  calculated. 

In  general,  age-test  urinalysis  strategics  trade  off  predictability  for  reduced  tail  area.  Here  we 
mean  tlie  tail  area  of  the  distribution  of  the  number  of  tests  in  a  fixed  time  period.  Age-test 
strategies  provide  fewer  people  not  tested  withui  1  year  and  fewer  people  tested  excessively  during 
1  year.  Age-test  strategies  are  also  more  predictable,  have  lower  variance  in  the  number  of  tests, 
and  as  a  result  are  subject  to  gaming  by  drug  users. 


Both  the  NRC  alternative  A  and  the  SCE  prc»ccss  have  some  undesirable  propcrtics.These 
strategics  involve  large  differences  in  the  testing  rates  between  people  tested  within  the  past  year 
and  those  who  were  not  tested.  This  implies  that  once  tested  there  is  a  high  probability  of  not  being 
tested  again  within  1  yeai.  These  probabilities  are  0.74  for  NRC  alternative  A  and  0.48  for  the  SCE 
process.The  SCE  process  is  such  that  almost  one  half  the  tests  every  month  arc  given  to  people  who 


Currently,  extensions  to  the  age-test  model  are  under  development.These  include  drug  usage 
patterns  and  their  impact  on  the  probability  of  detection  .This  will  help  quantify  the  trade-off 
between  predictability  and  the  reduced  tail  area  in  the  distribution  of  the  number  of  tests  mcntictted 
above. 
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Derivation  of  Optimal  Conditions 


This  appendix,  dcvclopis  tlic  solution  to  the  following  problem: 

1 


Minimize  n,  - 


'  1  - 1  +  (9/^2  -  ^d-l‘id)  '  Pc/*  i 


Subject  to  Jij  +  ,  -  ^1^3  —  1  nj  <  a 

Pd*  1 

/>,  >  p,  /  •»  1, 2, . . 


where  0  <  a,  P  <  1.  Clearly,  the  problem  is  equivalent  to 
Minimize  1  +  +  ■  •  •  ^d- 


91^2  •  ••  9d-  \9d 


Pa*i 


..  ^1^2- • -^dl^id  - 

subject  to  - - - ^1  < 

•'  Pd*  1 


P/>P;  /  - 1,2  J. 


(Al) 


(A2) 

(A3) 


For  given  value's  of  </a  and  pj  ^ the  objective  function  (A2)  is  laigest  at  -  9?  - ...  -  .  i 

H<moe,  we  wish  to  solve: 

Maximize  x-  — ~ — 

Pd*  1 


1-P 

(A4) 


Subject 


to - - -  <  Ot 

ii/(i-P/  ']i(i-(i-  p/'yp]+^ 


p./>P 


(A5) 


where  (A5)  is  the  result  of  manipulating  (A3)  and  noting  that  1  +(1  -  P)  + . . .  +  (1  -  P)*^' '  -  (1  -  (1-  p)^/ 
p.  Finall}',  solving  (A5)  for  a  we  have  the  equivalence 


Maximize  x 


(id 
Pd*  1 


Subject  to  x< _ r _ , _ 

-  (l-ajd-pT^ 


1-(1-P)‘' 

P 


y(a,  P) 


Pd>^- 


A- 1 


Choosing  jt  - ><a,  P)  with pj>fi  yields  a  solution.  This  reduces  to 


P<i*\ 


1 


Pd  _ 

><a,P) 


Pd>^ 


and  (p^,  p^^\)  must  lie  on  a  line.  Figure  A-1  illustrates  the  situation.  When  (a,  P)  -  (aj,  P,),  there 
are  multiple  solutions  lying  on  L,  to  the  right  of  the  vertical  line  -  p.  In  such  a  configuration  a 
unique  solution  can  always  be  obtained  by  choosing  the  intersection  of  L,  and  the  45°  line  L  :  p^ 
,.  The  case  (a,  P)  -  (02,  P2)  implies 1  <Pj.  In  this  configuration  a  unique  solution  can  be 
obtained  by  letting  pj  -  p  and  taking  p^  +  ,  on  Lj  at pj  =  p  or p^  + ,  =  (1  -  PV  y((x,  P). 


Figure  A*l.  Geometric  interpretation  of  optimal  region. 

For  any  optimal  solution,  the  minimum  value  of  n,  is  obtained  by  substituting  yfa,  p) 
for  qj!  Pj  in  (Al)  and  setting  -  ^2  “  —  “  1  “  1  -  p-  This  yields 


rc,(a,P) 


_ 1 _ 

l+(l.p)  +  ...  +  (l-P)«'->+(l.p)«ly(a,P) 

p(l-«) 

1-(I-P)‘' 


A-2 
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